Doping graphene with metal contacts 



G. Giovannetti^'^, P. A. Khomyakov^, G. Brocks^, V. M. Karpan^, J. van den Brink^'^, and P. J. Kelly^ 

^ Instituut-Lorentz for Theoretical Physics, Universiteit Leiden, 
P.O. Box 9506, 2300 RA Leiden, The Netherlands 
^Faculty of Science and Technology and MESA^ Institute for Nanotechnology, 
University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands. 
^Institute for Molecules and Materials, Radboud Universiteit, Nijmegen, The Netherlands. 

(Dated: July 23, 2008) 

Making devices with graphene necessarily involves making contacts with metals. We use density 
functional theory to study how graphene is doped by adsorption on metal substrates and find that 
weak bonding on Al, Ag, Cu, Au and Ft, while preserving its unique electronic structure, can still 
shift the Fermi level with respect to the conical point by ^ 0.5 eV. At equilibrium separations, the 
crossover from p-type to n-type doping occurs for a metal work function of 5.4 eV, a value much 
larger than the graphene work function of 4.5 eV. The numerical results for the Fermi level shift in 
graphene are described very well by a simple analytical model which characterizes the metal solely 
in terms of its work function, greatly extending their applicability. 
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Recent progress in depositing a single graphene sheet 
on an insulating substrate by micromechanical cleav- 
age enables electron transport experiments on this two- 
dimensional system [ll[2]. Such experiments demonstrate 
an exceptionally high electron mobility in graphene, 
quantization of the conductivity, and a zero-energy 
anomaly in the quantum Hall effect, in agreement with 
theoretical predictions [3l HI O [H [7]. The spectacular 
effects arise from graphene's unique electronic structure. 
Although it has a zero band gap and a vanishing density 
of states at the Fermi energy, graphene exhibits metallic 
behavior due to topological singularities at the iC-points 
in the Brillouin zone [3l [4] where the conduction and 
valence bands touch in conical (Dirac) points and the 
dispersion is essentially linear within ±1 eV of the Fermi 
energy. 

In a free-standing graphene layer the Fermi energy co- 
incides with the conical points but adsorption on metallic 
(or insulating) substrates can alter its electronic prop- 
erties significantly [SEHIOlIIIlIIlinilllllS]. Since 
electronic transport measurements through a graphene 
sheet require contacts to metal electrodes [U [121 HSl HZ] 
it is essential to have a full understanding of the physics 
of metal-graphene interfaces. In this paper we use first- 
principles calculations at the level of density functional 
theory (DFT) to study the adsorption of graphene on a 
series of metal substrates. The (111) surfaces of Al, Co, 
Ni, Cu, Pd, Ag, Pt and Au, covering a wide range of work 
functions and chemical bonding, form a suitable system 
for a systematic study. 

Our results show that these substrates can be divided 
into two classes. The characteristic electronic structure 
of graphene is significantly altered by chemisorption on 
Co, Ni and Pd but is preserved by weak adsorption on 
Al, Cu, Ag, Au and Pt. Even when the bonding is weak, 
however, the metal substrates cause the Fermi level to 
move away from the conical points in graphene, resulting 



in doping with either electrons or holes. The sign and 
amount of doping can be deduced from the difference of 
the metal and graphene work functions only when they 
are so far apart that there is no wave function overlap. At 
the equilibrium separation, the doping level is strongly 
affected by an interface potential step arising from the 
direct metal-graphene interaction. 

Based upon the DFT results, we develop a phenomeno- 
logical model to describe the doping of graphene, taking 
into account the metal-graphene interaction. The model 
uses only the work functions of graphene and of the clean 
metal surfaces as input. For a given metal substrate, it 
allows us to predict the Fermi level shift in graphene with 
respect to the conical points i.e., both the type and con- 
centration of the charge carriers. The model also predicts 
how metal work functions are modified by adsorption of 
graphene. 

Some details of how DFT ground state energies and op- 




FIG. 1: (Color online) The most stable configurations of 
graphene (a) on Cu, Ni and Co (111) with one carbon atom 
on top of a metal atom (A site), and the second carbon on a 
hollow site (C site) and (b) on Al, Au, Pd and Pt(lll) in a 
unit cell with 8 carbon atoms and 3 metal atoms per layer. 
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FIG. 2: (Color online) Band structures of graphene absorbed 
upon Al, Pt and Co (111) substrates. The Fermi level is at 
zero energy. The amount of carbon pz character is indicated 
by the blackness of the bands. The conical point corresponds 
to the crossing of predominantly pz bands at K. Top panels: 
graphene on Al and Pt; bottom panels: the majority and 
minority spin bands of graphene on Co. Note that on doubling 
the lattice vectors (for Al and Pt), the K point is folded down 
onto the K point of the smaller Brillouin zone. 



timized geometries are calculated for graphene on metal 
(111) surfaces are given in Ref. [l^. We fix the in- 
plane lattice constant of graphene to its optimized value 
a = 2.445 A and adapt the lattice constants of the met- 
als accordingly. The graphene honeycomb lattice then 
matches the triangular lattice of the metal (111) sur- 
faces in the unit cells shown in Fig. [l] The approxi- 
mation made by this procedure is reasonable, since the 
mismatch with the optimized metal lattice parameters 
is only 0.8-3.8%. We have verified explicitly that the 
structures shown in Fig. [l] represent the most stable con- 
figurations of graphene on the metal substrates studied. 
The equilibrium separations, binding energies and work 
functions are listed in Table [ij 



TABLE I: Calculated equilibrium separation deq of a graphene 
sheet from various metal (111) surfaces. The binding energy 
is the energy (per carbon atom) required to remove the 
graphene sheet from the metal surface. Wm and W are, re- 
spectively, the work functions calculated for the clean metal 
surfaces, and for free-standing and adsorbed graphene. 
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FIG. 3: (Color online) Calculated Fermi energy shift with 
respect to the conical point, AEf (dots), and change in the 
work function W — Wg (triangles) as a function of Wm — 
Wg, the difference between the clean metal and graphene 
work functions. The lower (black) and the upper (green/grey) 
results are for the equilibrium (^ 3.3 A) and a larger (5.0 A) 
separation of graphene and the metal surfaces, respectively. 
The solid and dashed lines follow from the model of Eq. ([T]) 
with Ac = for d = 5.0 A. The insets illustrate the position 
of the Fermi level with respect to the conical point. 



The results immediately show that the metals can be 
divided into two classes. Graphene is chemisorbed on Go, 
Ni and Pd(lll), leading to binding energies ~ 0.1 
eV/carbon atom and equilibrium separations (ieq ^ 2.3 
A. In contrast, adsorption on Al, Cu, Ag, Au and Pt(lll) 
leads to a weaker bonding, AE < 0.04 eV/carbon atom, 
and larger equilibrium separations, deq ~ 3.3 A. These 
results are in agreement with previous calculations and 
experimental data [8} [TO l [T2 l l2T| [22] . 

To identify the changes in the graphene electronic 
structure induced by adsorption, we calculate the band 
structures as illustrated in Fig. [2] for some typical exam- 
ples. When graphene is chemisorbed (on Go, Ni, and 
Pd) the graphene bands are strongly perturbed and ac- 
quire a mixed graphene- metal character. In particular, 
the characteristic conical points at K are destroyed, see 
the bottom panels of Fig. [2] When the interaction is 
weaker (Al, Gu, Ag, Au, Pt), the graphene bands, includ- 
ing their conical points at can still be clearly identi- 
fied; see the upper panels of Fig. [2] However, whereas 
in free-standing graphene the Fermi level coincides with 
the conical point, adsorption generally shifts the Fermi 
level. A shift upwards (downwards) means that electrons 
(holes) are donated by the metal substrate to graphene 
which becomes n-type (p-type) doped. 

For metal-graphene equilibrium separations, graphene 
is doped n-type on Al, Ag and Gu, and p-type on Au 
and Pt; the corresponding Fermi level shifts are plotted in 
Fig.[3j Because the work functions of graphene, Wq^ and 
of most metal surfaces, Wm^ differ, as soon as graphene 
interacts with a metal, electrons are transferred from one 
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FIG. 4: (Color online) Left: schematic illustration of the pa- 
rameters used in modeling the interface dipole and poten- 
tial step formation at the graphene-metal interface. Right: 
plane- averaged difference electron density An{z) = nM\G{z) — 
nyiiz) — nciz) showing the charge displacement upon forma- 
tion of the graphene-Pt(lll) interface. 



to the other to equilibrate the Fermi levels. A schematic 
representation is shown in Fig. [4] for the case of electron 
transfer from graphene to the metal. To a good approxi- 
mation, the graphene density of states (DOS) is described 
by D{E) = Do\E\, with Do = 0.09/(eV^ unit cell) for E 
within 1 eV of the conical points. Since this DOS is 
much lower than that of the metal, equilibrium is effec- 
tively achieved by moving the Fermi level in graphene and 
even a small electron transfer will shift the Fermi level 
significantly. A transfer of 0.01 electrons would lower the 
Fermi level by 0.47 eV. 

This electron transfer results in the formation of an in- 
terface dipole layer and an accompanying potential step 
AV. We can use the plane- aver aged electron densities 
n{z) to visualize the electron redistribution An{z) = 
^m|g(^) ~^m(^) — ^g(^) upon formation of the interface. 
As shown in Fig. |4j An{z) is localized at the interface. 
The sign and size of the interface dipole are consistent 
with the changes of the metal work function upon ad- 
sorption of graphene, see Table [l| 

Naively one would assume that graphene is doped 
with electrons if Wq > Wu and doped with holes if 
Wg < Wm' The crossover point from n- to p-type dop- 
ing would then be at I^m = ^g- The results obtained 
at the equilibrium separations of the graphene sheet and 
the metal surfaces {d ~ 3.3 A; see Fig. [S) show that this 
is clearly not the case. Instead, the crossover point lies at 
— Wg = 0.9 eV. Only when the graphene-metal sep- 
aration is increased significantly does the crossover point 
decrease to its expected value, as illustrated by the upper 
curve for d = 5.0 A in Fig. [s] This clearly demonstrates 
that the charge redistribution at the graphene-metal in- 
terface is not only the result of an electron transfer be- 
tween the metal and the graphene levels. There is also a 
contribution from a metal-graphene chemical interaction. 
Such an interaction, which has a significant repulsive con- 
tribution, has been found to play an important role in 




FIG. 5: (Color online) Fermi level shifts AEpid) as a function 
of the graphene-metal surface distance. The dots give the cal- 
culated DFT results, the solid lines give the results obtained 
from the model, Eg. ([!]) [25]. 



describing dipole formation when closed shell atoms and 
molecules are adsorbed on metal surfaces [23l |24] . 

The dependence of this interaction on the metal- 
graphene separation d is mapped out in Fig. [5] in terms 
of the dependence of the Fermi level shift AE-p on d. We 
use the parameters shown in Fig. [4] to construct a simple 
and general model with which to understand these re- 
sults. The work function of the graphene-covered metal 
is given by W{d) = Wu — AV{d) where AV is the po- 
tential change generated by the metal-graphene interac- 
tion. The Fermi level shift in graphene is modeled as 
AE-p{d) = W{d) — Wg- The key element is modeling the 
potential step AV = Atr(<i) + ^c{d) in terms of a "non- 
interacting" charge transfer contribution Atr driven by 
the difference in work functions and a contribution Ac 
resulting from the metal-graphene chemical interaction. 

The charge transfer contribution is modeled by a plane 
capacitor model as indicated in Fig.[4j Atj.{d) = aN{d)zd 
where a = e^je^A = 34.93 eV/A with A = 5.18 A^ 
the area of the graphene unit cell and N{d) is the num- 
ber of electrons (per unit cell) transferred from graphene 
to the metal (becoming negative if electrons are trans- 
ferred from the metal to graphene). Zd is the effective 
distance between the charge sheets on graphene and the 
metal. Zd < d diS most of the charge is located between 
the graphene layer and the metal surface as illustrated in 
Fig. [4j We model it as Zd = d — do with do a constant. 

Integrating the (linear) density of states of graphene 
yields a simple relation between N{d) and AE-p{d): N = 
:^DoAEp/2. Using the relations introduced in the pre- 
vious two paragraphs we can then express AE-^id) as 



AEF{d) = ± 



y^l^2aDo{d-do)\WM-WG-A,{d)\-l 



aDo{d — do) 



(1) 



where the sign of AEp is given by the sign of Wu — Wg - 



4 



Ac- The parameters do and Ac{d) turn out to depend 
only very weakly on the metal substrate. If we fit these 
quantities to the DFT results for one metal substrate, 
we can use them as universal parameters to predict the 
Fermi level shifts in graphene for all metal substrates. 
We use the DFT results obtained for graphene on Cu 
(111) to fix do and Ac{d), see Ref. [25]. Only the work 
function of the clean metal surface, I^m, and that of free- 
standing graphene, Wq^ are then needed to calculate the 
Fermi level shift. The accuracy of the model represented 
by Eq. ([T]) is demonstrated in Figs. [3]and[5l From AEf 
one can immediately obtain the work function W of the 
metal-graphene system, as well as the sign and concen- 
tration of the charge carriers in graphene, N. 

The critical metal work function Wu = Wo where 
the Fermi level is at the conical points of graphene, can 
be obtained from Eq. ([T]) for AEF{d) = 0. It gives 
Wo{d) = Wq + Ac{d). The contribution of the chemi- 
cal interaction term Ac depends strongly on the distance 
d between graphene and the metal surface. At a large 
distance d > 4.2 A, Ac < 1 eV and Wo{d) ^ Wq = 4.5 
eV, whereas at the equilibrium separation (ieq = 3.3 A, 
Ac ^ 0.9 eV and Wo{d) ^ 5.4 eV. This agrees with the 
DFT results shown in Fig. [3j The chemical interaction 
thus leads to a sizeable potential step at the equilibrium 
separation, which is downwards from metal to graphene 
as indicated in Fig. [4] The sign of this step and its in- 
sensitivity to the metal substrate are consistent with its 
interpretation in terms of an exchange repulsion between 
the electrons on graphene and the metal substrate [24 . 

In conclusion, we have used DFT calculations to study 
the doping of graphene induced by adsorption on metal 
surfaces and developed a simple model that takes into 
account the electron transfer between the metal and 
graphene levels driven by the work function difference, 
as well as the chemical interaction between graphene and 
the metal. The model extends the applicability of the 
detailed DFT results to the more complex systems en- 
countered in practical devices and suggests combinations 
of metal (strips) to be used to realize p-n junctions [26] , 
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